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The nonlinear rheology of a soft glassy material is captured by its constitutive relation, shear stress 
vs shear rate, which is most generally obtained by sweeping up or down the shear rate over a finite 
temporal window. For a huge amount of complex fluids, the up and down sweeps do not superimpose 
and define a rheological hysteresis loop. By means of extensive rheometry coupled to time-resolved 
velocimetry, we unravel the local scenario involved in rheological hysteresis for various types of 
well-studied soft materials. Building upon a systematic experimental protocol, we introduce two 
observables that quantify the hysteresis in macroscopic rheology and local velocimetry respectively, 
as a function of the sweep rate St^ 1 . Strikingly, both observables present a robust maximum with 
St, which defines a single material-dependent timescale that grows continuously from vanishingly 
small values in simple yield stress fluids to large values for strongly time-dependent materials. In 
line with recent theoretical arguments, these experimental results hint at a universal timescale-based 
framework for soft glassy materials, where inhomogeneous flows characterized by shear bands and/or 
wall slip play a central role. 



PACS numbers: 83. 60. La, 83.50.Ax, 83.50.Rp 

When submitted to an external stress, soft glassy ma- 
terials such as colloidal gels, clay suspensions, concen- 
trated emulsions, and foams, display a fascinating variety 
of behaviors because the applied strain may disrupt and 
rearrange the microstructure over a wide range of spa- 
tial and temporal scales leading to heterogeneous flow 
properties For more than a decade, flow dynamics 

have been probed by combining standard rheology, e.g. 
through the determination of the "constitutive relation" 
between the shear stress a and the shear rate 7, and local 
structural or velocity measurements 0, H[ • While much 
progress has been made on steady- state flow properties, 
the relevance of transient phenomena has been recog- 
nized only recently 0-0- Still, in practice, it can be 
argued that any experimental determination of the flow 
curve 17(7) is effectively transient since it is obtained by 
sweeping up or down 7 over a finite temporal window. 
In other words the measured flow curve coincides with 
the steady-state relation (7(7) only if the sweep rate is 
slow enough compared to any intrinsic timescale of the 
fluid. On the other hand, when the microstructure dy- 
namics are governed by long timescales, one expects hys- 
teresis loops in (7(7) measurements performed by sweep- 
ing up then down the shear rate (or vice versa). This 
phenomenon, referred to as "rheological hysteresis," has 
indeed been commonly observed in a host of complex 
fluids for about 70 years @, E| . However, to date, this 
ubiquitous signature of the interplay between timescales 
in complex fluids has not been quantitatively studied by 
means of local measurements. 

In this Letter, we use time-resolved velocimetry to un- 
veil the local scenario involved in rheological hysteresis 
in various types of well-studied soft materials. Build- 
ing upon a systematic experimental protocol, we intro- 



duce two observables, A a and A v , that quantify the am- 
plitude of the hysteresis phenomenon as a function of 
the sweep rate St^ 1 in macroscopic rheology and local 
velocity respectively. Both A a and A v go through a 
maximum with St, pointing to the existence of a char- 
acteristic timescale 8 for the microstructure dynamics. 
In thixotropic (laponite) suspensions and (carbon black) 
gels, 6 is of the order of several hundreds of seconds, 
while it becomes hardly measurable for simple yield 
stress fluids such as carbopol and concentrated emulsions. 
Velocity profiles allow us to understand this evolution 
by clearly differentiating a succession of homogeneous, 
shear-banded, and arrested flows depending on the fluid 
and on the sweep rate, thus providing a local interpreta- 
tion of rheological hysteresis. 

Experimental set-up and protocol.- Experiments are 
performed in a polished Plcxiglas Couettc geometry (typ- 
ical roughness 15 nm, height 28 mm, rotating inner 
cylinder of radius 24 mm, fixed outer cylinder of ra- 
dius 25 mm, gap e = 1 mm) equipped with a home- 
made lid to minimize evaporation. Rheological data are 
recorded with a stress-controlled rheometer (MCR 301, 
Anton Paar). Two flow curves are successively recorded, 
first by decreasing the shear rate 7 from high shear 
(7max = 10 3 s" 1 ) to low shear (7 m i n = 10~ 3 s" 1 ) through 
N = 90 successive logarithmically-spaced steps of dura- 
tion St each, and then by immediately increasing 7 back 
from 7 m i n up to the initial value 7 max following the same 
N steps in reverse order. In general the downward and 
upward flow curves, (7down(7) and <7 up (7) do not coincide 
and define a hysteresis loop [see Fig. [TJa) as an example] . 
Simultaneously to the flow curves, the azimuthal velocity 
v is measured as a function of the radial distance r to the 
rotor, at about 15 mm from the cell bottom, and with 
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a spatial resolution of 40 /im by means of ultrasonic ve- 
locimetry [Io|, 11 1. Velocity data are then averaged over 



the duration St of each shear-rate step at 7 and the cor- 
responding velocity profiles Vdown(7,?") and i> up (7,r) are 
normalized by the rotor velocity vq ~ je to allow for a 
direct comparison of flow properties at widely different 
shear rates [see insets in Fig. [If a)] fl2| . 

At 7 max the flow reaches a steady state within the 
time interval St, even for the fastest sweep rates, so that 
starting from high enough shear rates ensures a well- 
defined and reproducible initial condition. Our choice 
of 7min and N results from a compromise between a 
good sampling of the shear rate and a reasonable total 
duration. By monitoring the viscoclastic moduli prior 
to and after each experiment, we checked that evapora- 
tion and/or irreversible (chemical) aging of the sample 
were negligible even for the slowest sweeps (St = 300 s, 
2NSt = 15 h). Finally, approximating our step-like pro- 
tocol by a continuous sweep, the equivalent sweep rate is 
d(log7)/di = 1/nSt, where n = N/ log(7 max /7 min ) is the 
number of steps per decade. In the present protocol we 
fix n = 15 and identify the sweep rate with St -1 , while 
keeping in mind that nSt is the actual control parameter 
that probes the effect of structural and flow timescales 
in rheological hysteresis (see Fig. 1 in supplemental ma- 
terial). 

Laponite suspension.- Clay suspensions are well known 
to exhibit large rheological hysteresis [l]| . Here we focus 
on laponite samples prepared by mixing ultrapure water 
with 3% wt. of laponite powder (Rockwood, grade RD) 
and 0.3% wt. of hollow glass spheres (Sphericel, Potters) 
acting as acoustic contrast agents [10( . Figure Q] shows 
the hysteresis loops for St = 3, 10, and 100 s. To quantify 
their areas, we introduce the following observable: 
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where Aa(j) = o- up (j) — 0d O wn(7)- Note that the loga- 
rithmic sampling of 7 gives an equal weight to low and 
high shear rates. Strikingly, when the sweep rate is de- 
creased, i.e. when St is increased, A a goes through a 
maximum at St* ~ 25 s [Fig. [2(d)] . 

To uncover the local scenario underlying this non- 
monotonic behavior of A a , we turn to time-resolved ve- 
locity profiles recorded simultaneously to the shear-rate 
sweeps (see also movies in the supplemental material). 
For small values of St, the laponite suspension is quickly 
"quenched" from high shear rates to lower ones. Veloc- 
ity profiles are all linear both on the way down and on 
the way up even at the smallest shear rates [Fig. [TJa)]: 
the laponite suspension is not given enough time to re- 
structure and remains fluid throughout the whole cycle. 
The fact that the downward and upward flow curves do 
not superimpose is thus linked to small microstructural 
changes in the fluidized suspension (e.g. local fluctua- 



FIG. 1: (color online), (a)-(c) Flow curves a vs 7 of a 3% wt. 
laponite suspension obtained by first decreasing 7 from 10 3 
to lCP' 5 s _1 in 90 logarithmically-spaced steps of duration St 
each (black symbols), and then increasing 7 over the same 
range (red symbols). Each plot corresponds to a different 
time interval per step: (a) 5t = 3 s, (b) 10 s and (c) 100 s. 
Insets: velocity profiles inside the gap recorded at the same 
shear rate during the downward (black) and upward (red) 
sweeps. Velocity data are normalized by the rotor velocity vo 
and the vertical scale goes from to 1. The position r within 
the 1 mm gap of the Couette cell is measured from the rotor 
(r — 0). (d) Hysteresis loop area A a denned by Eq. ((TJ) vs St. 
(e) Area A v defined from the velocity profiles by Eq. ((2| vs 
St. Solid lines are lognormal fits of the A a and A v data. 



tions in the concentration or in the size of colloidal ag- 
gregates) responsible for differences in viscosity only |3[ . 

For intermediate values of St, the velocity profiles along 
the two flow curves now strongly differ [Fig.[IJb)]: while 
they remain linear during the downward sweep, they be- 
come inhomogeneous during the upward sweep and ex- 



3 



hibit shear banding over a large range of shear rates. 
Here, the suspension is given enough time to restructure 
leading to an arrested band subject to physical aging 
close to the fixed wall at r = e, that progressively disap- 
pears as the shear rate is increased. This contributes to 
enhance the size of the hysteresis cycle which thus mainly 
results from the competition between physical aging and 
shear rejuvenation [14|. 

Finally, for large values of St, flow arrest also occurs 
along the downward sweep [Fig. 0Jc)]: shear banding is 
observed for 7 ~ 0.1-1 s" 1 until the system experiences 
total slippage at the rotor (i.e. v = everywhere in the 
bulk) for 7 < 10~ 2 s _1 . This fully arrested state per- 
sists on the upward sweep up to much higher shear rates 
(7 ~ 10 s _1 ) and gives way to a homogeneously sheared 
state above a small interval of shear rates that corre- 
spond to decreasing shear stresses and to inhomogeneous 
flows. In the arrested state, the shear stress remains 
roughly constant, which tends to decrease the low shear 
rate contribution to the loop area, hence the decreasing 
behaviour of A a for St > 25 s. 

In order to rationalize the difference between the veloc- 
ity profiles recorded along the upward and downward flow 
curves, we consider Av(j,r) = 1^(7, r) - t>down(7, r), 
and integrate it over the whole gap and then over the 
whole range of explored shear rates in a way similar to 
Eq. (TTJ: 
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FIG. 2: (color online). Same as Fig. [TJfor a 1% wt. carbopol 
microgel. Solid lines in (c) and (d) are linear fits of the A a 
and A v data in semilogarithmic scales. 



As shown in Fig. [U(d), the evolution of this local ob- 
servable A v with St directly reflects that of the area A a 
extracted from the sole global rheology and its maximum 
is reached for a similar value of St* ~ 25 s. This confirms 
that the non-monotonic behavior of A a results from the 
flow properties at the mcsoscopic scale. This scenario is 
robust and does not depend significantly on the bound- 
ary conditions or on sample age since preparation (see 
Figs. 2 and 3 in supplemental material). We suggest that 
the characteristic time 9 = nSt* ~ 375 s is linked to the 
restructuration dynamics of the material. 

Carbopol microgel.- To further test the above idea, we 
turn to a simple yield stress fluid, namely a carbopol 
microgel, where restructuration is expected to be fast 
15l [iff . Results are reported in Fig. [2j The global ob- 



servable A a is of the same order as for the laponite sus- 
pension (the apparently smaller hysteresis in carbopol is 
due to the much larger vertical scale in the flow curves). 
Still, in the case of carbopol, A a is monotonically decreas- 
ing over the range of explored St [Fig. [2jc)]. Here again, 
the local observable A v follows the exact same trend 
[Fig. EJd)]. Velocity profiles reveal that the local sce- 
nario remains qualitatively the same whatever the sweep 
rate: homogeneous shear flow is observed along most 
of the downward sweep together with an ever-increasing 



amount of wall slip. For 7 < 0.1 s _1 the local shear rate 
vanishes and the flow becomes plug-like at roughly half 
the rotor velocity. On the way up, plug-like flow gives 
way to shear-banded profiles right after the stress maxi- 
mum and eventually to linear profiles at high shear rates 
[see insets in Fig. [Ha,b) and movies in the supplemental 
material], consistently with previous reports on startup 
experiments at constant shear rate 0,0]. However, as 
St is increased, the range of shear rates over which shear 
banding is observed gets narrower, leading to smaller val- 
ues of A v . In other words, the larger the time spent for 
each imposed shear rate, the closer the flow gets to its 
steady state, and thus the smaller the areas A a and A v . 

Discussion.- The above results hint at a unified pic- 
ture in which the behaviour of A& and A v is interpreted 
in terms of a single characteristic timescale 9 = nSt*: in 
laponite suspensions, we find 9 ~ 300 s while in carbopol 
microgels 9 is too short to be measured and only the de- 
creasing parts of A a and A v are observed. In Fig. [3] we 
explore the robustness of our findings by investigating 
other soft glassy materials. For three different concen- 
trations, thixotropic laponite suspensions [Fig. EJa)] and 
carbon black gels [Fig. 03h)] display a bell-shaped A a . 
The case of a concentrated emulsion [commercial may- 
onnaise, Fig. GHc)] appears as intermediate between the 
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FIG. 3: (color online). Area A a of the hysteresis loop as a function of the waiting time per point St for (a) laponite suspensions 
at 2.5% wt. (0), 3% wt. (□), and 3.5% wt. (o), (b) carbon black gels at 4% wt. (0), 6% wt. (□), and 8% wt. (o), and 
(c) a carbopol microgel at 1% wt (□) and a commercial mayonnaise (Casino) (o). The ■ symbols in (a) correspond to a 
3% wt. laponite suspension free of acoustic contrast agents and show that the addition of hollow glass spheres does not affect 
A a . In (b), carbon black particles and 1% wt. hollow glass spheres are suspended in a light mineral oil (Sigma, density 
0.838, viscosity 20 mPas), which allows experiments over 100 h (5t = 2 ■ 10 3 s) without any evaporation, (d) Normalized 
Aa-/A% data for all fluids except carbopol vs [St/ St*] 1 / A , where A%, St* , and A are the parameters of the lognormal fits, 
A a = A* a exp{ — [\og(St/5t*)/ A] 2 }, displayed as solid lines in (a)-(c). 



case of carbopol and the two previous materials: the loop 
area A a exhibits a maximum at the faster end of the ac- 
cessible range of sweep rates corresponding to St ~ 2.5 s 
i.e. 9 ~ 40 s. These bell-shaped curves are well captured 
by lognormal fits, leading to the master curve shown in 
Fig|3|d). This suggests that the characteristic time 9 con- 
tinuously grows when going from simple yield stress fluids 
(carbopol, emulsions) to highly time-dependent materials 
(laponite, carbon black), in line with recent theoretical 
arguments and numerical simulations [l5, 17 1. However 
the maximum value of A a and the width of the lognor- 
mal law have no straightforward dependence on material 
type or concentration and deserve further investigation. 
Whether or not 9 is linked to some intrinsic timescale gov- 
erning steady-state now properties also stands out as an 
open question. Still, interestingly, in both carbopol mi- 
crogels and laponite suspensions, inhomogencous shear- 
banded flows are observed concomitantly with sections 
of the upward flow curves where a presents a plateau 
[Fig. Hlb)] or decreases with increasing 7 [Fig. [3c) an< i 
[2ja,b)]. This is in agreement with a universal criterion 
proposed recently in Ref. [l8j which predicts transient 
shear banding to arise just after any stress overshoot dur- 
ing startup. 

To conclude our quantitative analysis of rheological 
hysteresis, based on well-defined protocols and observ- 
ables, unveils robust features in the transient flow curves 
that are linked to local flow properties such as shear 
banding and wall slip. As the waiting time per point 
St is increased, two regimes are observed. First, the ma- 
terial remains fluidized throughout the whole cycle and 
the area of the hysteresis loop increases with St contrary 
to the common belief that the slower the sweep rate, 
the smaller the hysteresis Q. Second, the area even- 
tually decreases as the material is given enough time to 



restructure. A well-defined material-dependent timescale 
9 separates these two regimes and apparently vanishes for 
simple yield stress fluids. Our study thus provides exper- 
imentalists with observables that quantify the distance 
to steady state in measurements of the constitutive rela- 
tion: whatever the soft glassy system, sweep rates larger 
than 9~ 1 should be used in order to ensure that effects of 
long-lived transients and inhomogeneous flows are mini- 
mized. Clearly, spatially-resolved models are needed to 
further explain and complete our experimental findings. 
Soft-glassy rheology approaches d, [l8| and local fluid- 
ity models [13, Ell appear as ideal candidates to check 
whether rheological hysteresis could be envisioned in a 
universal framework as suggested by the present work. 
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SUPPLEMENTAL MATERIAL 

Rheological hysteresis 
in soft glassy materials 

Supplementary movies 

Five movies are provided as supplemental material to 
show the full data sets corresponding to Figs. 1 and 2 in 
the main text. The files names indicate the system under 
investigation (laponite suspension or carbopol microgel) 
and the value of the waiting time per point St. Veloc- 
ity profiles are shown only for 7 < 250 s _1 due to the 
limitation in the ultrasonic pulse repetition frequency to 
20 kHz, which sets an upper bound on the measurable 
velocities [10| . At very small shear rates noise becomes 
too large for velocities to remain reliable. Indeed, the 
amount of velocity data recorded at each shear-rate step 
of fixed duration St increases with 7 [loj . Therefore, at 
small 7, the lower statistics result in significant noise in 
the velocity profiles. 

In the movies, the data corresponding to the downward 
(upward resp.) sweep are shown in red (blue resp.). The 
velocities profiles measured during the downward sweep 
are replotted with smaller symbols together with the up- 
ward sweep data, in order to allow for a direct comparison 
of the two velocity profiles obtained at the same shear 
rate. Black squares at r = and r = 1 mm indicate 
respectively the rotor velocity (vq ~ -ye) and the stator 
velocity (v = 0). 



Influence of the number of steps within a cycle 




laponite suspension for different values of n, the num- 
ber of measurement points per decade in the shear rate 
sweeps. The data all nicely collapse on a single curve 
when plotted as a function of nSt. In the experiments 
reported in the main text n is fixed arbitrarily to 15 
and only St is varied. The collapse shown in Fig. [ljb) 
indicates that the true control parameter is nSt. There- 
fore, when interpreting the time interval St* at which A a 
reaches a maximum, one has to keep in mind that the cor- 
responding physical timescale is rather 9 ~ nSt* ~ 500 s 
in the case of Fig. [1] 



Influence of the boundary conditions 

We checked for the influence of the boundary condi- 
tions on both obscrvablcs A a and A v . Figure [2] reports 
results obtained on two different laponite suspensions 
with smooth and rough boundary conditions. "Smooth" 
walls refer to polished Plexiglas, as used throughout the 
main text. The typical roughness of smooth walls is 
15 nm. "Rough" boundary conditions refer to a sand- 
blasted Plexiglas Couette cell (roughness of about 1 /xm) 
with the same dimensions as the smooth cell. We could 
not detect any systematic influence of the boundary con- 
ditions on A a and A v . In particular, the positions of the 
maxima St* in Fig. [2ja) and (b) are not affected by the 
wall roughness although the values of the hysteresis area 
A a can be up to 30% larger or smaller over the whole 
range of investigated 7 depending on the boundary con- 
ditions. 

In Fig. El^c), A v is seen to be about twice as large in 
the rough cell as in the smooth cell. A close inspection 
of the velocity profiles indicates that the material tends 
to stick at the stator in a rough cell so that fully arrested 
profiles (with v = everywhere) persist for longer times 
in the upward shear rate sweeps, leading to larger values 
of A v . The larger scatter in the A v data of Fig. [2](d) for 
rough boundary conditions also seems to arise from more 
erratic slip phenomena than in smooth conditions. Still 
the general evolution of the velocity profiles during shear- 
rate sweeps with rough walls remains otherwise very close 
to that observed with smooth walls. 



Influence of the natural aging process 



FIG. 1: (a) Area A a of the hysteresis loop as a function of the 
waiting time per point St for different numbers of measure- 
ment points per decade n = 3 (■), 8 (•), 15 (A), 30 (o), and 
50 (+). (b) Same data plotted as a function of nSt. The gray 
curve is a lognormal fit of the rescaled data. Data obtained 
on a 3.5% wt. laponite suspension in a smooth Couette cell. 

Figure [1] shows the curves A a vs St obtained on a 



Laponite suspensions are known to exhibit strong ag- 
ing properties due to slow chemical and physical pro- 
cesses [20(. We tested the influence of the sample age 
on rheological hysteresis by preparing a laponite suspen- 
sion that was left at rest for a time t w before part of 
the sample was loaded into the smooth Couette cell and 
submitted to scries of shear rate sweeps as described in 
the main text. Figure [3] confirms that aging processes 
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FIG. 2: (a,b) Area A a of the hysteresis loop as a function 
of the waiting time per point St in smooth (□) and rough 
(■) boundary conditions and (c,d) corresponding area A v de- 
duced from the velocity profiles for laponite suspensions at 
(a,c) 2.5% wt. and (b,d) 3.5% wt. Solid lines are lognormal 
fits. 
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FIG. 3: Area Aa- of the hysteresis loop as a function of the 
waiting time per point St in a 3.5% wt. laponite suspension 
for different ages t w since sample preparation: t w = 1 day (o), 
10 days (□), and 40 days (•). Experiments performed in a 
smooth Couette cell. The sample is free of acoustic contrast 
agents. The gray curve is a lognormal fit of the whole data 
set. 



